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Расчётная модель двухслойного пористого подшипника конечной длины 
с учётом анизотропии пористых слоёв и нелинейных факторов” 


А. Ч. Эркенов, М. А. Мукутадзе, В. С. Новгородова, Т. С. Черкасова 


Приводится асимптотическое решение по степеням относительного эксцентриситета задачи гидродинамиче- 
ского расчёта пористого двухслойного радиального подшипника конечной длины. Данное решение основы- 
вается на полных нелинейных уравнениях Навье — Стокса и на уравнении Дарси. Отличие предложенной 
расчётной модели состоит в том, что она учитывает анизотропию проницаемостей пористых слоёв, а также 
наличие источника смазки. Рассматривается случай, когда смазка принудительно попадает через поры 
внешнего пористого слоя. Найдено поле скоростей и давлений в смазочном и пористом слоях. В результате 
получены аналитические выражения для основных рабочих характеристик подшипника. Дана оценка влия- 
ния анизотропии проницаемостей пористых слоёв, нелинейных факторов, а также влияние источника смазки 
на основные рабочие характеристики подшипника. 

Ключевые слова: двухслойный подшипник, анизотропия проницаемостей, несущая способность, сила тре- 
ния, коэффициент нагруженности, коэффициент трения. 


Введение. Как известно [1, 2], устойчивый жидкостный режим трения в пористых подшипниках 
можно обеспечить только регулярной подачей смазки в зазор. Анализ существующих работ [1—5], 
посвящённых расчёту пористых подшипников, работающих под давлением питания, показывает, 
что в приведённых здесь расчётных моделях не учитывается анизотропия проницаемости пори- 
стых слоёв одновременно в радиальном и окружном направлениях. 

Для расчёта двухслойных пористых подшипников, работающих под давлением питания, 

необходимо представить коэффициент проницаемости в виде непрерывной функции, зависящей 
от радиальной и окружной координат. Однако учёт анизотропии только лишь в окружном направ- 
лении [3, 4] не позволяет представить коэффициент непроницаемости таким образом. Ниже нами 
приводится решение рассматриваемой задачи в нелинейной постановке при учете зависимости 
проницаемости пористых слоёв от радиальной и окружной координат. 
Постановка задачи. Рассматривается установившееся течение вязкой несжимаемой жидкости в 
зазоре двухслойного пористого радиального подшипника конечной длины. Подшипник с неодно- 
родным пористым слоем на его поверхности считается неподвижным, а шип вращается с угловой 
скоростью ©. Смазка в зазор пористого подшипника переменной проницаемости в осевом направ- 
лении подается под давлением питания. 

Поместим начало цилиндрической системы координат г, 9, 2 на оси подшипника на равном 
расстоянии /от его концов (1). 

Тогда уравнения контуров шипа и подшипника можно записать в виде 

с:Гг=Ь, с.:Г=Ь+В, с:Г=О+В, П=В,-В, 


(1) 


с. :г=а.(1+Н), Н = 56050 — Е? 5? 0, 


е 
где  — толщина пористого слоя; а — радиус шипа; Ь — радиус подшипника; ==-; е — 
а 


эксцентриситет. 


* Работа выполнена в рамках инициативной НИР. 
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Рис. 1. Схематическое изображение радиального подшипника с многослойным пористым вкладышем и пористым шипом 


Проницаемости пористых слоёв зададим таким образом, чтобы на границе раздела 
пористых слоёв они принимали одинаковые значения 


‚ , 

к = де“"®ы + (8), к-де” “+ (0). (2) 
Здесь А — заданная постоянная величина; безразмерные параметры Л, >0 и Л, >0 характери- 
зуют распределение проницаемостей пористых слоёв в радиальном направлении. Функцию К (6) 
по аналогии с законом изменения формы смазочной плёнки зададим в виде К =-А'с0$0, где 


А’/А = =" <1. Также предполагается, что =’и = — малые параметры одного порядка = =пП*='. 


Основные уравнения и граничные условия. Будем исходить из безразмерных уравнений 
движения вязкой несжимаемой жидкости в смазочном слое и в том числе в пористых слоях 
вкладыша, а также из уравнения неразрывности. 


2 2 2 2 2 
нь [91] “| 1 др бу 10, 10 бу у: 20% 
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и, щи, д 1 бр 0%, 19, 10%, 0%, 
Ве ч, + + = № + 
(1-а) 0925 0 ГОГ г“ 09 02 (3) 


би, 10 0, ц, 
+ + +—=0, 
ог г 00 02 г 
- 9Ф 10 т1оФ ох оКоФ 10 Ф 
К (г, 6). |-— + и + = 
0г° Гог г” 00° 02 дг 0г г 00 09 
- ФЕ 10Е 10?Е @Е] ОКЕ 10, дЕ 
К, (г,0)- * а + = 
г’ Гог г? 09° 02 ог ог г 09 69 








Й 





Здесь и,, ч,, у, — безразмерные компоненты вектора скорости; р — безразмерное гидродина- 


7! 


мическое давление в смазочном слое; Ф и Е — соответственно безразмерные гидродинами- 
ческие давления в пористых слоях; & и К, — безразмерные проницаемости пористых слоёв. 
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Размерные величины г,2,9,,ц,ч,,Р,Ф,Е и К,,К, связаны с безразмерными Г, 2, 0,, 
0,0, Р, ФЕ, К, К, соотношениями 
Г=Ш, 2=1, К=АК, К =АК, ц,=0а0,, и, =©аб, 
ы@ар + ф_ Нав Е _ Зав = (4) 


= ЕР: Ф, 
(р-а) Ь-а Ь-а 





В дальнейшем знак ^ у безразмерных переменных опускается. 

Система уравнений (3) решается при следующих граничных условиях. 

1. На контуре < выполняется непрерывность давлений [В = Ф) ‚ а компонента вектора 
скорости и, определяется законом Дарси. Остальные компоненты равны нулю. 

2. На границе раздела пористых слоёв 


фе. 
ог ог 
3. | „„=Рь, где р, — закон подачи смазки. 


4. На контуре с при 2 =+у давление равно р„, где р, = р, / р” : 
5. На поверхности шипа выполняются следующие условия: 





ц, (а+аН)=ч, | _ +[5=] аН +... = -=5т6, 
г=а аг ео 


(оон) -щ, .+[58) аН +...=1, (5) 
ии (а+аН) =и _ +5") аН+...=0. 
Гг=а ОГ ры 


Асимптотическое решение задачи. Установив закон подачи смазки как 
р, =с(2*-\?)+ У (а, с050 + Р, $т6)=”, (6) 
П=1 
решение системы (3), удовлетворяющее выше приведённым граничным условиям, будем искать в 
виде: 
р=А: (2? -\?)+Ра +Р(",6), ФЕМ. (2? -У)+ Ра + Е, (г,0), Е М. (2? -У)+ Ра +, (г,6), о 
у, =и(г,9), ч,=и(г, 0), ц, =и/(г,0)2 
Исходя из вида граничных условий функции Р, К, Ф, и(г,6), ч(г,9) и/(г,9) будем искать в 
виде рядов по степеням параметра = : 
Р=Р, (Г) +ЕР (^,0)+..., К, =Ф, (г) + =Ф, (п,0)+..., К, = Е, (г)+ЕР (п,9)+..., 
и=и (Г) + Е, (г,0)+..., у=ц, (г) + =, (г,0)+..., и=и, (г) + =, (г,0)+... 
Для определения коэффициентов разложений (8) с точностью до членов О(=^) придём к следую- 
щей системе уравнений и граничных условий к ним: 


(8) 








Ре| № 4 |_ 1 @ ‚9% 1% Ш 
За“ ЗЕ (1-а % 4 Гч г?" 
Ве тб 


а г | т га г’ 
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ещо ни ме |2 к : 














0 7-00 0 ро т’ 0 
ди, 10 и, -о, д°Ф, ‚10%, , 1 0°Ф, __^, 0, 
де г 00 г 9г? гг г? 06? г. бР 
2 2 
225 108 108 __№0Е. (10) 
07? гг г? 00? га 
Граничные условия запишем в виде: 
и ‚ (а) =0, ц, (а) =1, и ‚(а)=0, № (1)=А, № (1) =0, 
и (И =-ф(1-а)Ф, (1), (1) =0, и о а), $, (1)=Р, (0, 
Фо (В,) = А, (В,), Ф (В, = © (В), Фо (Вь ) = ра - Ра, М(Вь) = с (11) 
и, (а, 0) = -т6, и, (а, 9) =-и, (а о и 9) =-и,, (а)асоз0, 
и/ (1,0) =0, и, (11а) =-+(1- в) 52 +4 (1-а)Ф, (1]п' 0$0, 
=1 
ОФ, : 
и, (1,а) =-ф(1-а)—* 59! ‚ Р (1,0) =Ф, (1,0) =0,- (В,,0) =а, со$@ + В, те, 
Ф, (В,,9) =^ (8В,,0), Ф, (В,,9) = ^/(В,,0). (12) 
Г. Г. АБ мк 
Здесь В, =1+-№, В. =1+—, ф= „во 
| И (6-а} 


Исходя из вида граничных условий, решение системы (10) для первого приближения бу- 
дем искать в виде: 


и =и (г) со +, (г) тб, и, =ц,, (г) с0$ 9 + ци, (г) 6, 
и =и/ (г) со$9-+ии,, (г)эт@, Р =Р, (г)с0$0+Р, (г) те, 
и =и, (г) 05 9 + и, (г)зтб, Ф, =Ф,, (г)с0$0 +Ф,, (г) т6, 
Е =Е, (г) со50 +, (г) эт6. 
Подставляя выражение (13) в уравнения (10) и граничные условия (12), получим: 


" 1 , 2. 2 1 , В. ' , \ 2% 
и + и 2 и: 2 Ч? = 2 Р:: +ке 1 + Чои 1 + — Ир = Ч: ' 
Г Г Г (1 = а) Г 


(13) 
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2 1 | , : 4 2% 
И) +— Ч: — И? — И: м: 2 Р. +Ве ИоИ12 + Чо - — Ч: - Ч: |, 
Г Г (1-а) Ё Г 
2 1 1 у и у 
и +2 и = РР, +Ве | щи, +1, +4, +“. +- и 
1 11 Зы (а? ^ 12 091 И и 
ие 2 ей ь Тр +Ве| ии", + И! и + + 
о ао а а] Е 12 12 р ее 
" 1 , 1 , , 4 
И +-И - >И: = Ке Ио! + Иа + ИИ: |, 
Г Г Г 
" 1 ’ 1 ’ ’ 4 2 
И/ 2 И м = Ке | И‚и> + И, + И + АИ И/, |, 
1 и: 1 и 
и ++ =0, 4-2, +-В =0 
Г 8 а в , 
1 1 А 1 1 А 
Фи, + ро ‚Фи = - и. ‚ Ф, +-Ф,, - —Ф., = -ФЬ. (14) 
Г Г Г Г 
па) = 0-0 Е ч,, (а) = -, (а)а, 


1 
гы (9 0, и’, (а) = ии, (а)а, и’ (а) =0, 
ии (1) =-+(1-а)Ф, (1)- (1-а)ц, (1) п”, ц. (1) =-ф(1-а)Ф,, (1), 
ч,: (1) =-+(-а)Ф,, (1), ч, (1) =-(-а)Ф,, (1), 
и’. (1) =0, и’, (1)=0, р, (1)=0, р, (1)=0, 
Ф., (1)=0, Ф,, (1)=0, А: [В =а (В) =, 
ЕР: (В,) ыы Ф,, (В,), Е (В,) — Ф,, (В,), Ра (В,) = Ф., (В,), Е (В,) = Ф,, (В,) . (15) 
Заменяя в выражениях (9), (11), (14) и (15) производные слагаемые конечноразностными пред- 
ставлениями, получим систему алгебраических уравнений, которая решается методом Гаусса — 
Зейделя. 
Определение основных рабочих характеристик подшипника. Определив поле скоростей и 
давлений в смазочном слое, можно перейти к определению основных рабочих характеристик 
подшипника. Для составляющих вектора поддерживающей силы К, и К), а также для момента 
трения получим следующие выражения: 
=пура?ЬО/ =пура?6О/ 
о ИНВЕА Е „ (а), 8 Я ы (а), 


ЕЕ 
м си о) 


(16) 


т Ь а 


Основными рабочими характеристиками рассматриваемого подшипника являются: коэффициент 
нагруженности с, коэффициент сопротивления вращению &, коэффициент трения Ё Они опреде- 
ляются по формулам: 


м(1 г Е - 
= ‚ М= |2 +82, Е= =>. ТЯ 
- о А т ны, 1-а с чй 


Кроме того, параметрами рассматриваемого пористого подшипника, влияющими на его работо- 
способность, являются: 

— постоянная проницаемость стенки вкладыша, прилегающей к смазочному слою (характеризует- 
ся параметром 4); 
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— радиальный относительный эксцентриситет =" = Г. ; 
а 


— параметр п’ё, обусловленный переменной составляющей проницаемости в окружном направ- 


лении; 


ы ь +в 
— толщина вкладыша, характеризуемая безразмерной величиной В. = Г. : 





2. Й 
— отношение толщины пористых слоев г 
1 


| 
— длина подшипника / характеризуемая безразмерным параметром у = р р 


— параметры ро, р,, а,, ,, обусловленные наличием давления питания; 

— число Рейнольдса (Ве); 

— параметры Л, и ^Л»›, характеризующие распределение давления в радиальном направлении в 
пористых слоях. 





10° 10° 10° 10 | ‘и 


Рис. 2. Зависимость коэффициента нагруженности от параметра ц: Л, = Л, =О; а=0,998 ; р, -р, =0,5 ; В» = 1,6; 


1—у= 0,5; 2 — у = 1; сплошная линия — Вл = 1,1; пунктир — В1 = 1,2; пунктир с точкой — В! = 1,4 


27 


1-а 


10 





—3 —2 —1 
10 10 10 1 у 


Рис. 3. Зависимость коэффициента трения от параметра ф: Л, =Л, =0; а=0,998 ; р, -р, =0,5 ; В» = 1,6; 


1—у=2; 2 —у = 0,5; сплошная линия — Вл = 1,1; пунктир — В1 = 1,2; пунктир с точкой — В! = 1,4 
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Рис. 4. Зависимость коэффициента сопротивления от коэффициента проницаемости: Л, =Л, =0; а=0,998 ; 
р. -рР, =0,5 ; В = 1,6; 1—у=2;2 —у= 13 —у= 0,5; сплошная линия — В! = 1,1; пунктир — В+ = 1,2; 


пунктир с точкой — В; = 1,4 


Заключение. Результаты численного анализа, приведённые на рис. 2-4, показывают следую- 
щее. 

1. При значении параметра # ==, т. е. в случае, если проницаемость пористого слоя по 
окружности меняется по тому же закону, что и форма смазочной плёнки, подшипник по несущей 
способности обладает свойством подшипника двойного действия. 


2. При Л, =Л, =0, т. е. в случае однослойного однородного пористого вкладыша в зави- 
симости коэффициента нагруженности от параметра с при значениях у <1 наблюдается ярко вы- 
раженный максимум при ф=10`". 

3. При Л, =Л, =0 и у>1 с увеличением значения параметра ф значение коэффициента 
нагруженности уменьшается. Наиболее резкое снижение значения коэффициента нагруженности 
наблюдается при ф =103. 

4. При Л, =Л, =0, уе [1; 2] с увеличением значения параметра ф значение коэффици- 


ента трения увеличивается. При этом резкое увеличение коэффициента трения имеет место при 
значениях ц > 10-3. При значениях щ < 10-3 коэффициент трения практически остается постоян- 
ным с увеличением параметра ф. 

5. При Л, =Л, =О; у = 0,5; Ве [1,1; 1,2] в зависимости коэффициента трения от пара- 
метра ф наблюдается ярко выраженный минимум при ф = 1. 

6. Расчеты показывают, что при Л, * Л, , т. е. в случае двухслойного пористого подшип- 


л в т в 
ника, при значениях ги < и т. > 1 подшипник обладает более высокой несущей способностью, 
1 и 


К в 
чем при — >1и -2 <1. 
^, В, 
7. Приведённые выше (в пунктах 1—5) выводы остаются в силе и при Л, + Л», т. е. и в слу- 
чае пористого двухслойного подшипника при непрерывном изменении проницаемости пористых 
слоёв в окружном и радиальном направлениях. 
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Таким образом, на основании вышеизложенного можно сделать следующий вывод. При 
учёте анизотропии проницаемости пористых слоёв в окружном и радиальном направлениях мож- 
но обеспечить повышенную несущую способность подшипника при достаточно низком значении 
коэффициента трения. 
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ОЕЗТСМ МООБЕЁ ОЕ РОЧВЕЕ-ГАУЕК РОКОЦ$ ЕТМТТЕ-ГЕМСТН ВЕАВТМС \/ТТН 
ВЕСАВО ТО АМТ$ОТВОРУ ОЕ РОВОЦ$ ГАУЕВ$ АМО МОМЫМЕАВ ЕАСТОВ$” 


А. С. Егкепо\, М. А. МиКикад?е, У. $. Моудогодо\ма, Т. $. СпегКазоуа 


Тре азутрюйс 5о/ивоп ю те ргоМет оп те пуадгодупати/с сасшайоп оЁ {те рогоиз аоиЕ-Евуег габа! Етйе- 
епоЁёй Беаппд /5 ргомаеа оп {Те ессеттсйу гайо 1еуеб. Тре зо/ийоп 15 Базеа оп Ви! поп/теаг Мамег—юоКез 
едиайоп5, апа оп Рагсу5 едиаНоп. Тйе @5ИпсНоп оЕ {те зиддееа аея/дп тодеЕ! 15 {па # а!ои/з Гог пе регтеа- 
Бу агизоЁ гору о! 1е рогоиз [ауегз, ап ауайа Бу оЁ а дгеазтд 5оигсе. А сазе и/пеп дгеазта Юго деЁ 
{ргоидр рогез о! ап ехегпа! рогоис !ауег 15 сопз!Аегеа. А уеосйу апа ргеззиге Пе т ве 1ибтсапЕ апа рогоиз 
[ауегх 15 Юипа, Аз а гезш, апа/уйса! ехргез5опз Гог {пе тат регогтапсе (а о! {те Беаппд аге обеатеа. Тве 
еПесЕ оЁ те регтеа Бу апогору оЁ рогоиз !ауег5, поп/пеаг ГаСюг5, аз ИЕ! аз {Те тПиепсе о! а дгеазта 
50игсе оп {Те Бахс регогтапсе свагасетеНс$ оЕ те Беаппд аге ема/иае4. 

Кеуигога5: доиМе-вуег Беаппд, регтеа Шу аги/зоору, Беагпа сарасйу, Ейсвоп Гогсе, оайта соей«септе #/С- 
оп а&огЕ 


* Тне гееагс {5 аопе мт {Пе #тате оЁ {Не паерепдепе В&р. 
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